NYU-TH/00/03/08 
SO (10) Theory of R-parity and Neutrino Mass 



Charanjit S. Aulakh^, Borut Bajc^ 2,3 \ Alejandra Melfo*- 4 -*, Andrija Rasin*- 5 ^ and Goran 

Senjanovic^ 6,2 -' 
W Dept. of Physics, Panjab University, Chandigarh, India 
( 2 ) Department of Physics, New York University, New York, NY 10003, U.S.A. 

( 3 ) J. Stefan Institute, 1001 Ljubljana, Slovenia 
(4) (j AT, Facultad de Ciencias, Universidad de Los Andes, Merida, Venezuela 
^ Dept. of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599, U.S.A. 
^International Center for Theoretical Physics, Trieste, Italy 

Abstract 

We study the Higgs sector of a SO(10) grand unified theory which predicts 
exact conservation of R-parity at all scales and incorporates the see-saw mech- 
t^J- ' anism. We find possible intermediate scales and light states compatible with 

the constraints coming from the running of the gauge couplings. Such a pat- 
. tern could lower the SO(10) breaking scale, allowing the d = 6 proton decay 

operators to be comparable in magnitude to the d = 5 ones. 
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I. INTRODUCTION 



It is well-known that SO(10) grand unified theory offers an appealing framework for the 
unification of quarks and leptons and their forces. It is the minimal theory which unifies 
a family of fermions in a simple irreducible spinorial representation. The left-right (LR) 
symmetry in the form of charge conjugation is a finite gauge transformation and thus 
automatically built into the theory. Last but not least, it naturally incorporates the see-saw 
mechanism 0] which provides a rationale for the smallness of the neutrino mass. 

However, in its non-supersymmetric realization it fails to provide a "canonical" see-saw 
formula ||. It is here that supersymmetry (for a recent review on supersymmetric grand 
unified theories see for example 0) plays an important role and in certain cases the see- 
saw can take the canonical form. We study this issue in the context of the SO(10) theory 
and find, unfortunately, a situation similar to the ordinary case. More important, the see- 
saw mechanism actually gives us the low energy effective theory of supersymmetry. What 
happens is that R-parity remains an exact symmetry at all scales. This is a fundamental 
result which guarantees the stability of the lightest supersymmetric partner (LSP), an ideal 
candidate for the dark matter of the universe. This is true of any renormalizable theory of 
the see-saw mechanism based on the spontaneous breaking of B— L symmetry H. 
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Without unification one can not predict the value of the right-handed neutrino mass 
scale. In this sense SO(10) is an ideal theory especially in its supersymmetric version. It 
incorporates all the above features and it helps narrow down the range in which the right- 
handed neutrino mass lies. 

In this paper we construct a complete SO(10) theory of see-saw mechanism and R-parity. 
After a careful study of symmetry-breaking patterns, we compute the particle spectra and 
perform the unification analysis. We find a plethora of states (often carrying color or a 
large electromagnetic charge) whose masses, due to supersymmetry, could in principle lie 
much below the associated symmetry breaking scales. The reason for this is the violation 
of the survival principle, which we discussed at length in There we have coined the 
term "survival of the fittest" for this phenomenon. It is due to the absence of some quartic 
couplings in the potential (i.e., the absence of some trilinear couplings in the superpotential) , 
as is often the case in supersymmetric models. 

The rest of this paper is organized as follows: in the next section we discuss the salient 
features of the SO (10) theory with a renormalizable see-saw and the possible symmetry 
breaking patterns. In Section III we give detailed analysis of the particle spectrum and 
in Section IV we present the unification constraints. In Section V the phenomenological 
and cosmological consequences are discussed at some lengths. Section VI is devoted to the 
summary. The technical and computational details of symmetry breaking are left for the 
Appendices. 

II. THE SO(10) SUPERSYMMETRIC THEORY WITH THE 
RENORMALIZABLE SEE-SAW MECHANISM. 

Supersymmetric SO(10) models have been studied at length ||, but almost exclusively 
with the non-renormalizable version of the see-saw mechanism. More precisely, one chooses 
one (or more) pair of Higgs fields in the spinorial representation 16 and 16 whose VEVs 
induce B-L breaking and the mass for the right-handed neutrino through the d = 4 terms: 
m UR ~ (1616)/Mpj. The disadvantage of this program is that then R-parity is broken at a 
large scale Mr, and thus one needs additional, often ad- hoc symmetries to understand the 
smallness of R-parity breaking at low energies. 

Our motivation is orthogonal to this. We wish to have a theory of R-parity and this 
points immediately to the renormalizable version of the see-saw mechanism. To see this, 
recall first that under R-parity p —>■ p, p —>■ —p (where p stands for particles and p for 
sparticles of the supersymmetric standard model) can be written as (S is the spin) 

R= (-lf( B ~V+ 2S = M(-1) 2S , (1) 

where the so-called M parity is obviously equivalent to R (the factor (— 1) 2S becomes 1 for 
the physical Hamiltonians and only scalars with S = are allowed to have non-vanishing 
VEVs). 

One of the most appealing aspects of the supersymmetric SO(10) is that M is a finite 
gauge symmetry, since under M 

16 -16 , 10 10, (2) 
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and all other representations built out of the fundamental 10, such as 45, 54, 126, etc. are 
even. The symmetry in (fj) is simply C 2 , where C is the center of SO(10), so that under 
it 16— > i 16, 10— > — 10. This points strongly towards using a 126-dimensional Higgs for 
the breaking of B— L and the see-saw mechanism (for previous analysis of SO (10) models 
with 126 representations see [PH-I2]|). One drawback of these representations is their huge 



contribution to the SO (10) (3 function, so that the Landau pole is reached very soon above 
Mx- We will be interested only in the physics below or at the unification scale, where the 
couplings are perturbative. 

We wish to construct a renormalizable SO(10) theory with a see-saw, and this requires 
the minimum set of Higgs representations which break SO(10) down to the MSSM: 

S = 54 , A = 45 , E = 126 , E = T26 . (3) 

Although SO(10) is anomaly-free, as is well-known, one has to use both E and E in order 
to ensure the flatness of the D-piece of the potential at large scales 3> Mw- 

What are the possible channels of symmetry breaking and why do we need all the above 
fields? Of course, one possible channel strongly encouraged by the MSSM unification con- 
straints is through SU(5). Another possibility, much less studied, is the intermediate Pati- 
Salam (PS) and/or LR scale. For this purpose, it is useful to know the decomposition of 
the above fields under the PS gauge group SU(2) L x SU(2) jR x SU(4) C 

S = 54 = (1, 1, 1) + (1, 1, 20) + (3, 3, 1) + (2, 2, 6) , 

A = 45 = (1, 1, 15) + (3, 1, 1) + (1, 3, 1) + (2, 2, 6) , 

S = 126= (3, 1,10) + (1, 3, 10) + (2, 2, 15) + (1, 1, 6) , 

E = 126= (3, 1, 10) + (1, 3, TO) + (2, 2, 15) + (1, 1, 6) . (4) 

One may think that it is redundant to take both A and S together with E and E. After 
all either (1, 1, 1) or (1, 1, 15) could produce the first stage of symmetry breaking, down to 
the PS or LR symmetry, to be followed by (1,3,10) in E. Notice that the (1,3,10) field 
contains a color singlet component A c which couples to the right-handed neutrino v c in 16, 
and through (A c ) = Mr produces the see-saw mechanism. Since v c is an SU(5) singlet, so is 
obviously (A c ). Thus the possible (126) = (A c ) ^ can break SO(10) only down to SU(5). 

Clearly, either (A) ^ ^ (E) or (S) ^ ^ (E) suffices to break the SO(10) symmetry 
all the way down to the MSSM. So why is it that we need both A and S fields in the 
renormalizable version of the theory which we study here ? 

To see this, we need to write down the most general superpotential containing S, A, E 
and E 

W = ^Tr S 2 + ^Tr S 3 + ^Tr A 2 + ATr A 2 S 
2 3 2 

+m s EE + risY> 2 S + rf s Z 2 S + r^EEA . (5) 

The above form is only symbolic, details are given in Appendix A. 
It is trivial to see why A, E and E cannot suffice. From 

F A = m A A + rj A Y^ (6) 
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and the fact that (EE) preserves SU(5), it is clear that (A) must be SU(5) invariant, too. In 
order to achieve more breaking, one may resort to non-renormalizable operators of the type 
Tr A A /M, Tr (EE) 2 /M, etc. (where M is some new large scale). We prefer to include another 
field S] after all if (S) 3> (A), we will get the non-renormalizable terms after integrating out 
the heavy fields. 

It is also straightforward to see why S does not suffice together with E and E. The point 
is that the rjs and 7} s terms cannot produce any interaction between the VEV acquiring 
fields (1, 1, 1)5 and (1, 3, 10)^: there is no singlet in the tensor product 10 x 10 in SU(4); the 
same argument applies for E. Thus, although (S) can break SO(10) down to PS, the lack 
of interactions for the (1, 3, 10) and (1, 3, 10) fields forces their VEVs to vanish. 

In short, in the renormalizable theory one needs all the above fields. One can envision 
the following physically interesting patterns of symmetry breaking, which we label (a) and 
(b) 



SO(10) 



(S) 



M x 



SU(2) L x SU(2) H x SU(4) C 




SU(2) L x SU(2) R x SU(3) c x U(1) SL SU(2) L x U(l) fl x SU(4) C 



(a) 



M R \ <£> , <£) 



SU(2) L xU(l) y x SU(3) C 




(b) 



M w 



SU(3) c x U(l) t 



(7) 



It turns out that both of these chains are quite interesting in the sense that they lead to 
a plethora of "light" states, i.e., states whose masses lie below the corresponding symmetry 
breaking scale. 

In order to study any of the chains, we need the conditions for the F and D flatness at 
the scales M w . This again is discussed at length in Appendix A, here we address only 
the salient features. The F-flatness equations are 

1 2 

F(i,i,i)s = m sS + ^ss 2 + -A(a 2 - b 2 ) = , 



F (i,i,i5) A = m Aa + 2\as + -rj A aa = , 
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^(i,3,i)a = m Ab - SXbs + -riAVO = , 
^(i,3,io) s = (T [m s + ri A (3a + 26)] = , 

F (li3>10) _ = a [m E + r] A (3a + 26)] = , (8) 

where 

s = ((1,1, l)s) , a=((l,l,15) A > , 6=((1,3,1) A ) , 

<7=<(1,3,10) E ) , a= ((1,3,1%) • (9) 

All other fields have zero vacuum expectation values (see Appendix A). Notice that the 
choice of the two chains of symmetry breaking depends on the ratio of a and 6. In both 
cases it is assured that s is the largest VEV. 
Next, we can imagine two possibilities: 

(a) s>a>a = a>6~ a 2 / s , 

(b) s>6>a = o->a~ a 2 / s , 

which correspond precisely to chains (a) and (b) respectively (notice that a = a is necessary 
for D-flatness, see Appendix A). This can be achieved by paying the usual price of fine- 
tuning, an issue which is beyond the scope of this paper. In case (a) this implies 

2 

m j4 + 2As~ — <s, (10) 
a 

which then ensures 

2 

6~ — <cr (11) 

s 

(it is important to keep in mind that 6 can never vanish). In case (b), of course, the 
conditions are obtained interchanging the roles of a and 6, i.e. m A — 3As ~ a 2 /b s. 

A comment is in order. In the F flatness conditions (§) we ignore the fields in 16 and 10 
dimensional representations. This is justifiable for the Standard Model non-singlet fields, 
but not for v c in 16. However, as we already noticed, v c is coupled to the SU(5) singlet in 
E which gets a VEV a. The F-equations will then give 

F* = a(v c ) = (12) 

guaranteeing a vanishing VEV for v c . This is a general feature of theories with a renormal- 
izable see-saw mechanism, as shown in 0. We readdress this important issue below, when 
we discuss the fate of R-parity. 

A few words about the low-energy sector of the theory. Besides the above fields we 
also need the usual 16-dimensional representations, which give the standard model fermions 
and sfermions and representations containing the electroweak Higgs doublets. For the sake 
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of minimality we choose for the latter the 10-dimensional representation. We leave their 
number open in order to be as general as possible. It is well known that the minimal theory 
with just one 10 predicts wrong fermion masses and no CKM mixing if one restricts oneself 
to the tree level and/or renormalizable interactions. The higher dimensional operators can 
easily introduce small flavor mixings and correct the light quark mass ratios. Recently it 
has been emphasized that the radiative corrections in supersymmetric models can also do 
the job |nj. We present our results for both the minimal model and the model with two 
10-dimensional Higgs supermultiplets, since the latter certainly works. 



III. PARTICLE SPECTRUM. 

Since we wish to determine the scales of symmetry breaking, we need the precise particle 
spectrum of both light (~ Mw) and heavy (^> Mw) states. The point here is that in 
supersymmetric theories one can not rely on the survival (or extended) principle, since 
the lack of cubic terms in the superpotential for many representations naturally suppresses 
many particle masses 0. More precisely, if there is an intermediate scale Mj besides the 
unification scale M x , then one finds the effective quartic couplings suppressed by Mj/M x 
and there are a number of states whose masses become of the order of Mj/Mx instead of 
Mi. This of course has a strong impact on the unification predictions. This will be made 
manifest in the examples discussed below. 

Case (a) 

It is a simple exercise to show the well-known fact that all the states in 54 become 
super heavy. So do most of the states in 45, 126 and 126. The only states which do not 
pick up a mass of order s ~ M x in 126 and 126 are the SU(2) triplets, namely (3, 1, 10)e, 
(1, 3, 10)s, (3, 1, lO)^, and (1, 3, lO)^- As mentioned above we must fine-tune the mass of 
the field (1, 1, 15) A (m A + 2Xs ~ 0). 

Next, we switch-on a ~ Mc and a ~ Mr with a ^> a. As is clear from (|7|), a ~ Mc 
breaks PS down to LR; it is the color singlet component in (1, 1, 15)^ that gets a VEV. 
Due to the supersymmetric version of the Higgs mechanism, color triplets in (1, 1, 15)^ get a 
mass of order a, but the color octet remains in principle much lighter. Due to the absence of 
the trilinear term A 3 , it can only get a mass of order M c /M x (through the mixing with the 
color octet in (1, 1, 20)^) or M\jMc (through the fine-tuning condition + 2As ~ a 2 /a). 

Regarding the fields in £ and S, notice first that the neutral component in (1, 1, 15)^ 
is a parity (charge conjugation) odd super-Higgs. In other words, the left-handed triplets 
(3, 1, 10)e and (3, 1, 10)^ get masses m s — 3ar]A ~ M c and decouple for lower energies. 
Similarly, the color triplet and sextet fields in their right-handed counterparts (1,3, 10)s 
and (1,3, 10) ^ get masses ~ Mc- 

Except for the doubly-charges states 5+ + and <5+ + , the rest of the fields in the color 
singlet components for the above representations get masses ~ Mr through the super-Higgs 
mechanism. The 5+ + and 5+ + fields, much as in LR supersymmetric models |14]], pick up 
their mass only through b ~ a 2 /s. 

The masses are summarized in Table 1 below. 
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State 


Mass 


all of S in 54 

all of A in 45, except (1, 1, 15)^ 

all of £ in 126 + S in 126, except SU(A) C decuplets 


~ M x 


(3,1,T0) S + (3,1,10)^ _ 
color triplets and sextets of (1, 3, 10)s and (1, 3, 10)^ 
color triplets of (1, 1, 15)a 


~ M c 


from the color singlets of (1, 3, 10) s and (1, 3, 10W 

O \ 7 7 J \ 7 7 / 2.1 


~ Mr 


color octet and singlet of (1, 1, 15)^ 


~ M x = Max [g, 


from the color singlets of (1, 3, 10) s and (1, 3, 10)^ 


~ M 2 = M 2 R /M X 



Table 1: Mass spectrum for the symmetry breaking chain SO(10) -?SU(2)^xSU(2)/jX 
SU(4) C ^ SU(2) L x SU(2) i? x U(l) B _ L x SU(3) C ^ SU(2) L x U(l) y x SU(3) C . 

Notice the 'light' states in the two last rows of the type we discussed above as the violation 
of the survival principle. Their small masses are the product of the lack of renormalizable 
interactions and fine-tuning conditions. 

On top of this there are of course three generations of fermions and sfermions which lie at 
M w or below. The right-handed neutrino supermultiplet is at M R , but it does not affect the 
running. In the minimal model the 10-dimensional Higgs multiplet splits into a light (~ My/ 
or supersymmetry breaking scale Mg) bidoublet, which provides the two MSSM doublets, 
and the superheavy (~ M x ) color triplet and antitriplet. If there are two 10s, then besides 
the two light doublets and all the superheavy colored states we will have two more doublet 
superfields with masses >> M w . The range for the heavy doublet masses depend on one's 
scenario for generation of the weak mixing angles. If these are to arise at the tree level from 
bi-doublet mixing induced by the vev of the (1,3,1) submultiplet of the 45, then one needs 
more than one 10-plet and it is necessary to assume Mh ~ b (where Mh is the coefficient 
of the 10 2 term in the superpotential), while the coupling A# ~ b/M x (where A# is the 
coefficient of the 54 10 2 term in the superpotential), in which case the extra doublets have 
mass ~ M R /M X . Else the mixing is negligible and the doublets are superheavy so that the 
theory effectively reduces to the case with a single 10. One may obtain realistic mixing even 



7 



with a single 10-plet if radiative corrections due to soft terms are appropriate 
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Case (b) Let us now turn to case (b) in (j^). This implies interchanging the values 
of the 45 VEVs a and h. The analysis proceeds along the same lines, so we just present the 
particle spectrum in Table 2. 



State 


Mass 


all of S in 54 

all of A in 45, except (3, 1, 1)^ + (1, 3, 1)^ 

all of £ in 126 + X in 126, except SU{A)c decuplets 


~ M x 


(3,1,T0) E + (3,1,10) E 
{!,"_ , 10)s and (1, , 1(% 
ujf from (1, 3, 1)a 


~ Mr 


color triolets and sinfdets from 

VVlVl UllUlVjllJ < 1 11V.1 0111C,1VjUO 11 Vlll 

(l,+,10) s and (1,-1% 


~ M c 


(3,1,1)a 


~ Mi = Max ^ ^ 


color sextets from 

(l,+,10) s and (1,-1% 


- M 2 - M 'f /Mx 



Table 2: Mass spectrum for the symmetry breaking chain SO(10) ^SU(2)xxSU(2)/jX 

SU(4) C ^ SU(2) L x U(l) fl x SU(4) C ^ SU(2) L x U(l) y x SU(3) C . The states in (1,3,10) 
and (1, 3, 10) were decomposed according to their T 3R number, for example (1, +, 10) denotes 
the component of (1, 3, 10) with T 3R = +1, etc. 

Notice again the survival of the fittest principle. These are the color sextets and left- 
handed triplet states. Again their small masses are due to the lack of renormalizable inter- 
actions or the fine-tuning condition [vtla — 3As ~ a 2 /b ~ Mq/Mr). 

In this case the heavy bi-doublets will have masses of order b ~ M R provided M H , X H 
are again chosen small. Else one may again use a single 10-plet and radiative effects. 



S 



IV. UNIFICATION OF GAUGE COUPLINGS 



Armed with the complete spectrum of the physical states we are now ready to perform the 
analysis of the unification constraints. Notice that this is a precise analysis without any ad 
hoc assumptions made regarding the spectrum at intermediate scales. At this stage we wish 
to have a rough qualitative estimate of the new mass scales and thus it is only appropriate 
to perform this at the one-loop level. In this way we check whether the assumptions (a) and 
(b) are consistent with the low-energy values of the coupling constants. 

We define the running coefficients as usual, 



1 



1 



bi , E 2 
— In — 

2tt E x 



(13) 



ai{Ex) a t (E 2 ) 

where i = 1, 2, 3 stands for the normalized hypercharge, SXJ(2) L and SU(3)c gauge couplings, 
respectively. 

For the multiscale case, which we are interested in, eqs. (0) become for the chain (a) 

— I Oi ' m— — h o- 'in-rrr- + °i ' iri 

atu 



2tt 



Oi(M z ) 



M s {2) M 2 ( . . 

W z + *i ln W s + b > ln W 2 + 



•(4) 



m Mi + 1 



M, 



In- 



M 



x - 



R 



Mi 



c 



(14) 



where Mi = Max[Ml/M c , M%/M x ], M 2 = M|/M x ; for case (b) one has to interchange 
Mr and Mc in flli]) as well as in the definitions of Mi and M 2 . 

Let us study the cases (a) and (b) defined in (0) separately, 
(a) SO(10)^SU(2) i xSU(2) i? x SU(4) c ^SU(2) L x SU(2) fi x U(l) B _ L x SU(3) C A -^ 
SU(2) L x U(l) y x SU(3) C . The coefficients are 



Energy range 








M z < E <M S 


-(41 + n)/10 


(19-n)/6 


7 


M s < E < M 2 


-33/5 


-1 


3 


M 2 < E < Mi 


-(57 + 3e)/5 




3 


Mi < E < M R 


-(57 + 3e)/5 







M R < E < M r 


-(45 + 3e)/5 







M c < E < M x 


-(191 + 3e)/5 


-(41 + 6) 


-34 



(15) 



where e counts the number of heavy Higgs doublets that kicks in at M R /Mx and (n + 1) 
is the total number of Higgs doublets below the supersymmetry breaking scale Ms- The 
detailed equations that follow from (]TJ]) specific to these models, as functions of e and n, 
are given in Appendix B. They can be solved to obtain the three mass scales Mx, Mc and 
Mr in terms of the unification coupling cxu- Consistency will determine Mi, which is the 
maximum of M R /Mc and M c /Mx- Notice that we have assumed in ( |I~5|) that Mi < Mr, 
which is not necessarily true. However equation ( |i~5|) shows that only b\ is sensitive to the 
scale Mr (all the other coefficients remain constant), while only 6 3 can feel the other scale 
Mi. For this reason it is irrelevant for the form of the renormalization group equations which 
of the two scales is smaller. 
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We have performed the calculation for values of the coupling constants a\[Mz) = 
0.01695, a 2 (M z ) = 0.03382 and a 3 (M z ) = 0.119 ± 0.003. 

For the case with just one 10 (e = 0), the unification scale M x turns out to be indepen- 
dent of a v . Mi is Mx/Mc in this case, and it is found, for n = 

1Og (^)^ 16 - 20±0 - 07 -°- 141Og (^) 

log (^) = 14.43 ± 0.07 - 0.25 log + 0.07^ 

log(J^|) = 13.71 ±0.12 -0.18 log +0.10^ (16) 

The lowest values are obtained for larger a 3 . With M s ~ 1 TeV and a\j ~ 1/5, we get 

log {^) - 14-77 ± 0.07 log ~ 14.21 ± 0.12 (17) 

Adding an extra 10-dimensional Higgs (e = 1) spoils the independence of M x - In this 

case Mi = M c /M x . We get, for n = 

log (^) = 14.38 ± 0.25 + 0.14 log (^) + 0.09^ 

log (^) = 12.77 ± 0.25 + 0.03 log (^) + O.lSc^ 1 

log(J^) =10.25 ±0.46 + 0.34 log + 0.28^ (18) 

We have the interesting result of M x being lower for lower intermediate scales M R . 
The smallest possible value of (a^/) -1 is dictated by proton decay constraints, requiring 
\og(M x / GeV) > 15.5. Again, raising the supersymmetry breaking scale has the effect of 
lowering the intermediate scales, and for Ms of order 1 TeV 



l0g (H) 



14.63 ±0.17 log (7^7) ^ 13.62 ± 0.29 (19) 



The lowest possible value corresponds to a 3 = 0.122, and gives otjj 1 = 9.5. 

This low value of Mr is very interesting for neutrino physics, and so we concentrate on 
this case, e = 1. Figure 1 shows that the three scales M x , Mq, Mr increase as a function 
of l/ajj, for M s = 1 TeV. It is interesting to see that the requirement of proton stability 
(M x > 10 15 - 5 GeV) guarantees the perturbativity of the unified coupling, l/a v > 12. 

In Figure 2 we plot our findings for M R and M c as functions of M s for a fixed value of 
M x = 10 15 ' 5 GeV. Notice that Mr is almost insensitive to Ms (and completely so for n — 1), 
unlike Mc- Of course, the value of Ms — 1 TeV is physically most interesting. Notice that 
there is some sensitivity to the precise value of 0:3. 

(b) SO(10) J ^SU(2) L xSU(2) i? x SU(4) C ^ SU(2) L x U(l)*xSU(4) c ^ SU(2) L x U(l) y x 
SU(3) C . 

The b coefficients are now (here e is the number of bidoublets from 10 with mass Mr) 
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Energy range 


b[ k > 






M z < E <M S 


-(41 + n)/10 


(19-n)/6 


7 


71 ,T 1 — 1 71 T 

M s < E < M 2 


-33/5 


-1 


3 


M 2 < E <M l 


-97/5 


-1 


-2 


M l <E <M C 


-97/5 


-3 


-2 


M c < E < M R 


-81/5 


-3 





M R < E < M x 


-(189 + 3e)/5 


-(43 + e) 


-30 



(20) 



It can be seen that the extra bidoublets do not play a significant role in this case, 
contributing only at scales above Mr. In the minimal case of e = 0, n = 0, we get Mi = 

Mlj Mr, and 

log (^) = 14.88 ± 0.11 - 0.13 log + O.OGc^ 1 

log (^) = 13.19 ± 0.11 - 0.241og (J^L) + 0.12a,} 1 
log(J^) = 13.04 ± 0.14 -0.171og(^)+0.13a^ (21) 

Again, lower intermediate scales are obtained for lower Mx- With Mx ~ 10 15 ' 5 GeV,in 
this smaller value of Ms is preferred. With Ms ~ Mz, 

log ( J^) ^ 14.47 ± 0.12 log ~ 14.35 ± 0.12 (22) 

the lowest value obtained for smaller a s ,and a\j is now 8.7. The case with the extra bidoublet 
gives essentially the same results. 

What we observed here is that there can be a separation of all the scales involved, 
although sometimes not very big. The presence of extra light states (in addition to those 
of MSSM) below M x , especially those that have non-trivial hypercharge, tended to spoil 
the succesful unification for too big separation of scales. However, as shown in (|H|), it is 
perfectly possible to push scales to the low intermediate range. Finally we remind that the 
results presented here are only a rough guide, since we have not included the two-loop and 
threshold effects. 



V. PHENOMENOLOGICAL AND COSMOLOGICAL CONSEQUENCES 

We now turn to the discussion of the resulting physics of the theory. The most important 
features are the exact R-parity at low energies, i.e. the MSSM as an effective low-energy 
theory and the "quasi-canonical" see-saw form with a rather light neutrino. 

i) R-parity. Here the situation parallels the one already found by some of us in the 
context of LR symmetric theories \TM. We showed that at the large scale (16) = (u c ) = 



and so R-parity is conserved. Next, as in any theory based on the renormalizable see-saw 
mechanism O, one can show that (z>) = too. This is valid at all energies all the way to 
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the electroweak scale and below (of course one has to assume (u) = at high scale, as with 
other fields that break the SM symmetry). 

It is worth repeating this simple, but important argument. It is well-known that in 
the MSSM we cannot break R-parity spontaneously. The problem is this would imply the 
existence of the Majoron J, since this also breaks lepton number spontaneously ||15|| , and 
the Z boson would then have an extra decay into a Majoron and its real partner. 

Now, after we integrate the heavy fields out, we are left effectively with the MSSM: at 
low energies all the effects of the new scale Mr must go as I /Mr. The same is true of the 
would-have-been Majoron J: its mass must be suppressed by I /Mr and Z can still decay 
into J + partner. Again, as in the MSSM, this is ruled out and we learn that R-parity can 
never be broken. This is an extremely important prediction, since then the LSP is absolutely 
stable. 

ii) Neutrino Mass. One of the main motivations behind the SO(10) unification is the 
natural implementation of the see-saw mechanism. We chose 126 and 126 Higgs multiplets 
in order to have a renormalizable see-saw which then gave us R-parity to be exact all the 
way down to the MSSM. 

The first thing we notice is that for both patterns of symmetry breaking the allowed value 
for Mr is pretty high (usually about 10 14 GeV or higher, although it can go down to 10 12 
GeV, as in ([TE|)). Second, supersymmetry helps us not to have uncontrollable VEVs from 
the left-handed triplet VEVs, and indeed the supersymmetry breaking effects are negligible 
for such high M R 

Further, as we emphasized in ||, the exact form of the see-saw is rather model dependent 
in super symmetric SO(10) theories. The canonical form for the see-saw for our model is 
spoiled by effective non-renormalizable terms in the superpotential of the form 

' A$ 2 A C (23) 



M 



x 



(where A and A c are left and right-handed triplets respectively, from 126 and $ is a bidou- 
blet from 10) which arise once the heavy 54 fields get integrated out. Such terms then 
produce nonvanishing VEVs for the left-handed triplets since they get generated from the 
term above with m^AA to give 

M x mz (A c ) 

where e = (A c ) 2 /Mx^s can be anywhere between 1 and 10~ 4 in the above chains of 
symmetry breaking. 

Thus the formula for neutrino masses is 

m^ife-hl)^- (25) 
m UR 

where / is the coupling of triplets and right-handed neutrinos. 

We see now that the e term can not be neglected compared to the canonical term coming 
from the neutrino Dirac masses, at least for the lighter generations. This might be a welcome 



addition to the otherwise rather constrained scenarios ||16|| . 
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Of course, in some specific models one can have this new non-canonical term in the 
see-saw formula vanishing. For example, this is what happens if the Dimopoulos-Wilczek 
mechanism [O is employed to solve the doublet-triplet splitting problem. Let us briefly 



repeat the argument |J. As mentioned above, the form (|23|) is obtained from the terms (H 
is the 10-dimensional Higgs) 



XsSJ: 2 + X H SH 2 + ^S 2 (26) 

after integrating out the heavy 5, for which the presence of the term X^SH 2 is crucial. How- 
ever, to solve the doublet-triplet problem this same term must be absent, which is obtained 
for example with an additional Z3 symmetry fl8fl . The Dimopoulos-Wilczek mechanism is 
then implemented by the new terms (among others) 

X W SH' 2 + XhwHA'H' , (27) 

where H and H 1 are 10's and A' a new 45 with the vev < A' >= diag(a', a', a', 0, 0) x T2. 

Since the XhSH 2 term is forbidden, only the SU(2)i doublets from H (but not from H') 
have zero mass. So the Dimopoulos-Wilczek mechanism at the same time solves the doublet- 
triplet splitting problem and gives the canonical see-saw. The price to pay is however the 
ad-hoc addition of a discrete symmetry and the inclusion of more extra fields. 

iii) Proton decay and related issues. The usual lore is that the unification constraints 
favor the minimal supersymmetric SU(5) grand unified theory. One gets the single scale GUT 
breaking Mx — 2 x 10 16 GeV, which makes the usual d = 6 contribution to proton decay 
quite small: 

r(p -> 7r°e + ) > 10 36±L V- (28) 



and the dominant decay mode p — > K + u comes from the d = 5 operators |I9|. It is also 
known that in SO(10) models, due to misalignment of the Higgs triplet and doublet Yukawa 
couplings and depending on the details of the underlying flavor theory, other decay modes 
to kaons may be dominant, such as p — > K°l + [SO . 



However, what we see in the types of models studied here, is that unification through 
alternate symmetry breaking channels is quite possible. Moreover even the "smoking gun" 
signature of supersymmetric models with kaon production in proton decay, may not be 
dominant. Namely, we saw that, contrary to the common belief, lowering of the intermediate 
scale Mr can cause lowering of the grand unified scale (at least in the case with more than 
one 10 Higgs multiplet). This then opens the possibility that the gauge boson mediated 
proton decay may be the dominant one and thus we could even have the possibility where 
in a supersymmetric model the dominant decay products be pions instead of kaons []! 

Let us also comment on other higher dimensional operators. We know that the dangerous 
d = 4 proton decay is absent, since R-parity is exact even at E < My/- Of course, R-parity 



x An important, although not dominant contribution of the p — > 7r°e + mode has been found in 
[21, but there it was produced by d = 5 operators. 
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cannot guarantee the absence of higher dimensional operators which break the baryon and 
lepton number. Since 16 4 is invariant under SO(10), one can have 



AW effective = C®^ + ... (29) 

Even if M = Mpi, unless C is extremely small, the above interaction would be danger- 
ously large. Of course, there is no proof that it must be present, but still, it would be nice 
to have some explanation for its absence. 

If one dislikes ad-hoc discrete symmetries, as we do (in spite of the euphemism that they 
can come from superstrings) , one natural possibility is to go to E(G) GUT. In E(6) the basic 
representation is 27 = 1650(10) + 10so(io) + 1 .50(10) > an d the dangerous interaction (27) 4 is 
not allowed. Of course, one must then split naturally the Higgs in 1050(10) from the new 
states D and D c (down-quark like) in 10. The study of this is beyond the scope of this work 
and is related to the construction of the E(6) grand unified theory. 

Of course, even without going beyond SO(10) or including nonrenormalizable operators, 
one gets proton decay mediating d = 5 operators ( P9"| ) as usual after integrating out the heavy 
color triplets and antitriplets, i.e. (1, 1, 6) from 10 (the same contribution as in SU(5)), but 
also from 126 |2l[ . 



VI. SUMMARY 

There are a number of well-known reasons that make SO(10) a popular grand unified 
theory. The two main ones are the grouping of families of fermions in 16-dimensional spinor 
representations, and the natural incorporation of the see-saw mechanism. The supersymmet- 
ric extension of the theory provides yet another important rationale, namely, the otherwise 
ad hoc discrete R-parity symmetry of the MSSM becomes a finite SO(10) gauge transfor- 
mation. As we showed in this paper, the symmetry remains exact throughout all the stages 
of symmetry breaking down to the electroweak scale. Thus, the low energy effective theory 
is completely determined: it is the MSSM with absolutely stable LSP. This result is a gen- 
eral property of the supersymmetric see-saw mechanism and a spontaneously broken B-L 
symmetry 0. 

An important phenomenological prediction of the theory regards the proton lifetime. 
Remarkably enough, the theory allows for an intermediate see-saw scale at the expense of 
lowering the unification scale compared to the single-step breaking of supersymmetric GUTs. 
As argued in ||, the existence of intermediate scales is due to the violation of the so-called 
survival principle, i.e. to the fact that there is a number of supermultiplets whose masses 
lie below the corresponding scales of symmetry breaking. The typical value of their masses 
is ~ M\jMx- We find that the lowest possible values of Mr are of order 10 13_14 GeV, and 
with Mx ~ 10 16 GeV, these states turn out to be out of experimental reach. 

Contrary to the usual belief, the low unification scale allows for the interesting possibility 
of proton decay being dominated by the usual d = 6 mode characteristic of ordinary GUTs, 
p — > e + 7t°. Now, this is tied up, as we said, to an intermediate see-saw scale in the range 
10 13_14 GeV. This is encouraging for the neutrino mass expectations. However, we must 
admit that (unlike the often present claims in the literature), one cannot really predict 
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neutrino masses. This is clear from the central formula (p5|), which shows that besides the 
usual canonical term proportional to the squared Dirac Yukawa coupling, there is another 
piece which cannot be predicted by the SO(10) symmetry itself. Clearly, there is enough 
freedom to accomodate the solar and atmospheric neutrino data, but no way of making 
honest predictions. 
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APPENDIX A 

We will denote SO(10) indices by i, j, .., SO(6) indices by a, b, and SO(4) indices by 
a, (3, ... With this convention, the field S in the 54 representation is represented by the 
symmetric, traceless second order tensor S^; the field in the 45 representation, A, by the 
antisymmetric second order A^-; and the fields S and £ in the 126 and 126 representations 
by the fifth-order antisymmetric tensors Syw m and Tiijkim, with self- and anti-self-duality 
conditions 

y . ]_ p y y V 

^-'ijklm ^tijklmopqrs*- J opqrs j '—'ijklm ^^ijklmopqrs^-'opqrs 

The superpotential is 

W = ^TrS 2 + ^TrS 3 + ^TrA 2 + \TrA 2 S 
2 3 2 

+ m s SS + + r] s Z 2 S + r^SEA , (30) 

where 



^^-'ijklm^-»ijklpA m p , SS ^Jijklm'^-'ijklm j S »S • 



The F-terms are then 

(•^Oji = m A^-ij + ^(A-uSij + SnAij) + -^■[^klmpj'^'klmpi ~ ^klmpi^klmpj} 

(Fsh = m s S l3 + \ S [(S% - ^TrS 2 ] + A[(A%. - ^TrA 2 ] 

_i_ ^[^_y y I VS y y 

~r ^klmpi^klmpj ~r ^ z -'klmpi z -'klmpj j 

( ^ 'y^) ijklm TTlS^-iijklm A VA[Aip^ipjklm A ~^^ijklmopqrsA t^Jtpq r s\ 
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1 



% 

ijklm ^T^ijklm VA[^-ip^-'pjklm, ^ijklmopqrs^ot ^-*tpqrs\ 

'^VS^^ip^-'pjklm ~~^£ijklmopqrsSot^tpqTs\ • (31) 

The VEVs for S and A that ensure the required symmetry breaking SO(10) — > 
SU(2) L xSU(2) fl xSU(4) c are 

S = sdiag(l,l, 1,-3/2, -3/2) x 1 2 , A = diag(a, a, a, b, b) x r 2 . (32) 

The VEV of £ and £ are a bit more complicated. In the basis where the Cartan 
generators are diagonal, the SU(2) £ xU(l)y xSU(3)c singlets in £ and £ get VEVs: 

(^13579) = > (^24681o) = ® ■ (33) 

One then goes to the basis in which (j32j ) is valid, to find the 2 5 components of £ and £ that 
get a nonvanishing VEV. They have the first index equal to 1 or 2; the second to 3 or 4, 
etc. Each have a VEV 2~ 5 / 2 a, with a factor of —i for each even index in £ and a factor of 
+i for each even index in £. 

The VEV directions are chosen to get the required symmetry breaking, and of course are 
not the only possible ones. To make sure that the vacuum chosen is not connected to the 
undesired ones by flat directions, one has to examine the mass spectrum: the flat directions 
will be associated with massless excitations. It is clear from the analysis in Section III that 
no such states exist for our vacuum. 

The VEV equations are then written as 

a[m a + r] A (3a + 2b)} = a[m a + r] A (3a + 26)] = (34) 

(m A + 2\s)a + ^aa = (35) 
(m A - 3\s)b + yffff = (36) 

4ms - ys] + \\{a 2 - b 2 ) = (37) 

Then if one fine-tunes (m A + 2Xs)a ~ M R , the first pattern (a) of symmetry breaking is 
obtained 

50(10) SU(2) L x SU(2) R x SU(4) C ^ SU(2) L x SU(2) R x SU(3) C x U(1) BL 

^ SU(2) L x SU(3) C x U(1) Y (38) 

with 

Mx = s ~ ms ~ m A ; Mc = a ~ mj ; M\ = aa (39) 
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and b ~ M^/M x . 

The second pattern (b) is obtained by instead fine-tuning (m^ — 3Xs)a ~ Mj|. Then 

50(10) ^ SU(2) L x SU(2) R x SU(4) C Q 5C/(2) X x l/(l) fl x SU(4) C 

^ SU(2) L x SU(3) C x C/(l) y (40) 

= s ~ ms ~ 3> M| = a(T > Mc = b ~ (41) 

and a ~ M C /M X - 

The contribution of the representations S(54), A(45), E(126) and E to the D-terms is as 
follows: 

Dij = (—ig u )(2S^ i Sj] k - 2A* k[i Aj] k + 5S^ mn[i S i]Mmn + 5E Mmn ^Ej] k i mn + ...) . (42) 

The contributions of the two index representations S and A may be written as commu- 
tators 

= {-i gu ){2[S\ S] + 2[A\ A] + , (43) 

but then eq. ( |3"2| ) for the vevs of the matrices S and A immediately give zero for their 
contribution to the D terms. 



Thus, it remains to see what the 126, 126 representations give. It is most convenient 
to work in the complex basis where the Cartan generators are diagonal (indices A,B = 
1,2. ..10) which is obtained from the usual Cartesian basis = 1,2. ..10) by the unitary 

transformation U = I5 X U2, where XJ% = -j= y | 1 A by Va = U T AiVi- 

It is easy to see that in this basis the 126 and 126 part of Dab = UAiUsjDij is 

Dab = {Sigu) (P*gcdef9g[a^b]cdef + ^gcdef9g[a^b]cdef) = D l AB — D l BA , (44) 
where Qab is the metric in the complex basis g — I5 x n and 

D\ B = (—5igu) (pGCDEFdGA^BCDEF — ^gcdef9gbEacdef) (45) 

can have non- vanishing vev (when the vevs of E, E take the form of eq. (|33|)) only if A,B 
form an even-odd conjugate pair of indices. Thus, without loss of generality, we have for A 
even and B odd 

< Dab >= -tgu(5\)g A B(W\ 2 - \a\ 2 ) + ... . (46) 

It follows that provided the magnitudes of a, o are equal, the D terms vanish at the first 
two stages of symmetry breaking specified by the vevs in eqs. (0), (|33|). 

The absence of flat directions characterized by holomorphic invariants formed from the 
fields involved in the high scale symmetry breaking (S, A, E, E) alone then follows post-hoc 
from the calculated spectra of these fields which do not include any massless modes (aside 
from those eaten by massive gauge bosons). The usual argument f2^JT4J] serves to invoke 



the protection of the soft terms against charge and color breaking flat directions involving 
the light field vevs. 
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APPENDIX B 



We now give in detail the renormalization group equations for each symmetry-breaking 
pattern. 



Case (a) M 2 = M 2 R /M X . 

For Mi there are two possibilities: 



1) M x = M 2 R /M C 



Then (Jj) becomes 



2n 271 ^.W i(2)\o , fouP) 1.(5) ^ 



Oi(M z ) au 

(~bf + bf + bf - b?)C + (-&?> + fcf > + , (47) 

where Y = ln(My /Mz) for any scale My. Using (|i5|) we get 



2, 2, /25-, ls+ 36 ± 6 ifl+ 146 c ,_215 ± 6 !x 



ai(M z ) a,7 V 10 
2tt 



a 2 (M z ) a v V 6 / 

27T 27T 

(AS + 6i? + 31C - 34X) . (50) 



a 3 (M z ) au 



2) Mi = Mg/M 



x 



In this case ( [Hj) becomes 
2tt 2tt 



ai(Mz) au 

(26 (3) _ 26 (4) + 6 (B, _ 6 j., )c7 + + ,(4) + 6f))x] ^ (51) 



and with numbers 



2. 2. ^5-« 5+ 36 + 6 £iJ+ 146 c _215 + 6 £x 



ai(M z ) ac/ V 10 

= -- ( 2 -^S + (2e)R + AOC - (41 + 2e)x) , (53) 
a 2 {M z ) au V 6 / 

= — - (AS + AOC - 37X) . (54) 
a 3 (M z ) au 



Case (b) M 2 = M 2 /M x . 



18 



We get the equivalent equations as (|14]), ( f47|) and (|5l[) just interchanging i? and C. We have 
again two cases: 

1) M 1 = M 2 C /M R 

2tt 2tt /25-n 112 _ 108 + 3e n 253 + 3e \ 

-o H C H ri A I , (55) 



ai(M z ) a n V 10 5 5 5 

2 " ^ /2 ^ n S + 4C+(38 + e)i?-(43 + e )X) , (56) 



a 2 (M z ) olu \ 6 
2tt 2vr 



a 3 (M z ) aj/ 



(4S + 8C + 30i? - 35X) . (57) 



2) Mi = M 2 R /M X 



2tt 



ai(M z ) 

2tt 
a 2 (M z ) 

2tt 
a 3 (M z ) 



2tt 



at/ 
2tt 

«t/ 



at/ 

2tt /25-n 



25 -n_ 112 108 + 3e o 253 + 3e 



6 



5 + (44 + e)i2 - (45 + e)X 



- {AS + 8C + 30R - 35X) 



(58) 
(59) 
(60) 
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FIGURES 
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FIG. 2. Mc and Mr vs. M s for Chain A, when Mx is fixed at its lowest allowed value 
(log(Mx /GeV) = 15.5) and e = 1, for n = and 1. Dotted lines are the same plots for lowest and 
highest values of 03 
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